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Physics is said to be universal when it emerges regardless of microscopic details. A prominent 
example is the Efimov effect, i.e., emergence of an infinite tower of three-body bound states obeying 
discrete scale invariance. Because of its universality and peculiarity, the Efimov effect has been 
subject to extensive research in chemical, atomic, and nuclear physics for decades. Here we show that 
collective excitations in Heisenberg magnets (magnons) also exhibit the Efimov effect by tuning an 
easy-axis exchange or single-ion anisotropy. We locate an anisotropy-induced two-magnon resonance, 
compute binding energies of three magnons, and find that they fit into the universal scaling law. 
We propose several approaches to experimentally realize the Efimov effect in quantum magnets 
and the emergent Efimov states of magnons can be observed with previously-used spectroscopic 
measurements. Our study thus opens up new avenues for universal few-body physics in solid state 
systems. 
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INTRODUCTION 

Sometimes we observe that completely different sys- 
tems exhibit the same physics. Such physics is said to 
be universal and its most famous example is the criti- 
cal phenomena pQ. In the vicinity of second-order phase 
transitions where the correlation length diverges, micro- 
scopic details become unimportant and the critical phe- 
nomena are characterized by only a few ingredients; di- 
mensionality, interaction range, and symmetry of the or- 
der parameter. Accordingly, fluids and magnets exhibit 
the same critical exponents. The universality in critical 
phenomena has been one of the central themes in con- 
densed matter physics. 

Similarly, we can also observe universal physics in the 
vicinity of scattering resonances where the s-wave scat- 
tering length diverges. Here low-energy physics is char- 
acterized solely by the s-wave scattering length and does 
not depend on other microscopic details. One of the most 
prominent phenomena in such universal systems is the 
Efimov effect, i.e., emergence of an infinite tower of three- 
body bound states obeying discrete scale invariance: 

_> A" 2 (n -> oo) (1) 

with the universal scale factor A ps 22.6944 [2]. Because 
of its universality and peculiarity, the Efimov effect has 
been subject to extensive research in chemical, atomic, 
and nuclear physics for decades after the first prediction 
in 1970 [2l|4]. In particular, the recent experimental re- 
alization in ultracold atoms greatly stimulated this re- 
search area [5HZ]- 

In spite of such active research, the Efimov effect 
has not attracted much attention in condensed matter 
physics so far. However, since the Efimov effect is uni- 
versal, it should emerge also in solid state systems. In 
this article, we show that collective excitations in Heisen- 
berg magnets (magnons) exhibit the Efimov effect by 



tuning an easy-axis exchange or single-ion anisotropy. 
We will locate an anisotropy-induced two-magnon res- 
onance, compute binding energies of three magnons, and 
find that they fit into the universal scaling law. We will 
also propose several approaches to experimentally realize 
the Efimov effect in quantum magnets. Because multiple- 
magnon bound states have been observed in various com- 
pounds with different methods [BU50"] , the observation of 
the Efimov states of magnons is within reach. Our study 
thus opens up new avenues for universal few-body physics 
in solid state systems. Also, in addition to magnetic re- 
alizations of Bosc-Einstein condensation [3~If[3~5] , the Efi- 
mov effect provides a novel connection between ultracold 
atoms and magnetic systems. 

TWO-MAGNON PROBLEM 

To demonstrate the Efimov effect in quantum mag- 
nets, we consider a Heisenberg spin model on a sim- 
ple cubic lattice with exchange (J z ) and single-ion (D) 
anisotropies: 

H = — — ^^(J S r+S + J z S r S r+ ^) 

r e 

-Dj2(S Z r) 2 ~Bj2 SZ r- ( 2 ) 
r r 

Here spin operators S^r = S^iiS^ and obey the usual 
commutation relations 

[S+, S",] = 2S*S ry , [5* S±] = ±S±5 r ,r> (3) 

and the identity 

(S+j 2s+1 = (5-) 2s+1 -0 (4) 

for a spin s representation. J, J z > corresponds to 
a ferromagnetic coupling and J, J z < corresponds to 



2 



an antiferromagnetic coupling. In the latter case, by 
rotating spins by ir along the z-axis, which leads to 
S^r — > — S 1 ^, only for sites with odd- valued x + y + z, 
we can choose J > 0. Therefore, J > is assumed from 
now on and is always understood as a sum over six 
unit vectors; Efc±i,±y,±£- 

A fully polarized state with all spins pointing down- 
ward, i.e., S$\0) = -s|0) and S~\0) = for all r, be- 
comes the ground state for sufficiently large B < 0. N- 
magnon excitations on this ground state are described by 
a wave function ^(ri, . . . , rjy) = (0| Jli^i I*)' which 
is symmetric under any exchange of coordinates (Bosc 
statistics) and satisfies the Schrodinger equation: 



EW(n,...,r N ) 
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*(n,...,r w ). (5) 



Here V| is an operator to displace to r ; + e and a 
constant energy shift —(6s J — 6s J z — 2sD + D + B)N 
is omitted in the right hand side. It is clear that J z > 
(< 0) acts as a nearest-neighbor attraction (repulsion) 
and D > (< 0) acts as an on-site attraction (repulsion) 
between a pair of magnons. Note that the effective N- 
magnon hardcore repulsion ^(r, . . . , r) = exists for s = 
(N - l)/2 according to {S-) 2s+1 = 0Q 

We now look for a two-magnon solution (N = 2) with 
a center-of-mass momentum K: 



*(ri,r 2 ) 



(6) 



where il = (ri + v%)f2 and p = ri — r 2 are center-of- 
mass and relative coordinates, respectively. The Bose 
statistics of magnons implies iPk(p) = ^k(~p) for the 
relative wave function. By treating the first term in the 
right hand side of Eq. ^ as a free part (H ) and the rest 
as an interaction part (V), the two-magnon Schrodinger 
equation can be brought into the Lippmann-Schwinger 



1 By setting r\ = 
N(N-l)D]V(r 



r N in Eq. jHJ, we obtain [E-3N(N- 1) J+ 
r) = for s = (TV — l)/2. Therefore, unless 
E - 3N(N - 1)J + JV(AT - 1)D = 0, the constraint . . . ,r) = 
is automatically satisfied. In particular, D plays no role for 
s = 1/2 due to *(r,r) = 0. 



equation: 

V>k(p) = 0K"(p) + (p 

= 0X(P) 



1 



r/a 



E-H a + i0+ 
dk cos(fc • p) 
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^{Jcos(f • e) - J z cos(fe • e)}ip K (e) - 2Dip K (0) 

. e 

(7) 



Here a is a lattice spacing, 4>k{p) is a solution to 
E4>k{p) = Hq4>k{p), and Ejc(k) is the energy of two 
non-interacting magnons given by 



(fc) = ^ sJ [2- 2cos( 



e) cos(fc • e) 



(8) 



By setting p = x, y, z and p = in Eq. Q, we ob- 
tain four coupled equations which determine the four 
unknown constants; tpKifi) and -0^(0), which in turn 
determine iPk(p) through Eq. Q. The same result in 
the isotropic case of J z = J and D = was obtained by 
the pioneering work in Refs. [3l)rt38| . 

A scattering resonance between two magnons is defined 
by the divergence of an s-wave scattering length a s where 
a two-magnon bound state appears from the threshold. 
a s can be extracted from the asymptotic behavior of the 
wave function at zero energy E — and zero center-of- 
mass momentum K = where 4>o{p) = C is a constant: 



0(|p|- 3 ). (9) 



lim tl) {p)\ =C\l-r—, 



By matching this asymptotic behavior with Eq. Q, the 
analytic expression of a s is obtained as 



b)] 



2-i + ^(i-^r)+3W[i-^-^(i-^)]' 



(10) 



where W ee ^ r(£) r( £) r(£) r(|l) « 0.505462 
results from one of the Watson's triple integrals [39] . 
Accordingly, we find that the two-magnon resonance 
a s — > oo takes place at 
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2.93654 



(11) 



for s = 1/2. Critical anisotropies for other spins are 
shown in Fig.[T] In particular, the two-magnon resonance 
for s = 1 is induced by the exchange anisotropy at 
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(D = 0), 



while it is induced by the single-ion anisotropy at 
D 12 

7 ~ swTi 



4.76874 



(Jz = J) 
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FIG. 1. (Color online) Critical anisotropy for the two-magnon 
resonance a s — > oo for s = 1 (solid), s = 3/2 (dashed), and 
s = 2 (dotted curve). There exist (s) no, one, or two s-wave 
bound state(s) at K — in the lower left, middle, or upper 
right region, respectively, separated by the two critical curves. 
An appearance of two degenerate d-wave bound states is not 
indicated here. 



or 

D 36T4 7 

-1 = 7vuF~i « 5 - 12703 ( J * = - J ) ( 13b ) 



TABLE I. Few lowest binding energies of three magnons in 
the s-wave channel at K = 0. Shown are results for s = 1/2 
and s = 1 with D — right at the two-magnon resonances 
induced by the exchange anisotropy [see Eqs. (Ill and ( 12 1]. 
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-2.00 x 10~ 6 


22.7 




CO 




22.6944 


TABLE II. Same as Table [I] but for s = 1 with J z = ±J 
right at the two-magnon resonances induced by the single-ion 
anisotropy [see Eqs. (|13a[) and (|13b[|]. 
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-8.88 x 10~ 6 


22.2 



for a ferromagnetic or antiferromagnetic coupling, respec- 
tively. Above these critical points, two magnons form a 
bound state at K = 0. Its binding energy E < is 
obtained by solving Eq. (J7| and determines the three- 
magnon threshold in Figs. [2] and [3] below. 



THREE-MAGNON PROBLEM 

With an anisotropy-induced two-magnon resonance 
being located, we now turn to a bound state problem 
of three magnons (N = 3). A bound state solution to 
the three-magnon Schrodinger equation can be obtained 
in a similar way to the previous two-magnon problem. By 
treating the first term in the right hand side of Eq. ^ as 
a free part (i?o) and the rest as an interaction part (V), 
the Lippmann-Schwinger equation for E < is written 
as 



*(ri,r 2 ,r 3 ) = (r 1 ,r 2 ,r 3 \ 



1 



E-H, 



■V\V). (14) 



Then we introduce a new parametrization of the wave 
function: 



x *(ri,r 2 ,r 3 ), 



(15) 



which describes three magnons with a center-of-mass 
momentum K in which two of them are separated by 



a distance p and the third one has a momentum k. 
The Bose statistics of magnons implies xi<{p',k) = 
Xk{— P',k). After a straightforward calculation, the 



Lippmann-Schwinger equation (14) can be brought into 

ri a 6k x cos( 2kl+ ^- K ■ p) 

_„ /a (27r/a) 3 E-E^kuK-k^-ksM) 



xk{p; k z ) 



^ {Jcos( ^ . e) - J z cos( 2fcl +^-^ • e)} XK (e;k 3 ) 



2^2{Jcos{ 



J z cos( kl+2 ^- K ■<§)} 



x fei) - 2D XK (0; fc 3 ) - 4D X k{0] fei) 



(16) 



where 



Ea(ki,h2,k 3 )= E J2 s J t 1 _ cos ( fc « ' «)1 ( 17 ) 



i=l,2,3 e 



is the energy of three non-interacting magnons. Finally, 



by setting p = x, y, z and p = in Eq. (16 1, we obtain 
four coupled integral equations which determine the al- 
lowed binding energy E < and the four unknown func- 
tions; XK"(ejfe) an d %fc(0;fc), which in turn determine 
XK(p;k) through Eq. fl6| . 

Because our purpose here is to demonstrate the Efi- 
mov effect of magnons, we from now on focus on the 
s-wave channel at K = where the Efimov effect is sup- 
posed to emerge. Accordingly, we look for bound state 
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FIG. 2. (Color online) Two lowest binding energies of three 
magnons for s = 1/2 near the two-magnon resonance (111. 
Dots are numerical results and the solid curve is from the 



universal theory in Eq. ( 18 1. The upper colored region repre- 



sents the three-magnon continuum and the cross (x) marks 
the critical point where a two-magnon bound state appears. 
For n — 0, non-universal regions where |a s |/a < 10 are indi- 
cated by the left and right gray areas. 



Xo(y; k z , k x , k y 

■ h h 



solutions where (i) Xo(x; k 

x i ky , k 

Xo(z] ky, k zi k x ), (ii) xo{&> k x , k y , k z ) = xo{x' : k x , k zi ky), 
(iii) Xo(0; k x , k y , k z ) is symmetric under any exchange 
among {k Xl k y , k z }, and (iv) all of them are even func- 
tions of k v [y — x, y, z). The resulting two coupled inte- 
gral equations with three variables ranging from to tt 
are solved numerically by discretizing each variable with 
the Gaussian quadrature rule. 

Tables |T] and [XT] show the few lowest binding energies 
of three magnons right at the two-magnon resonances 
located in Eqs. (|lT|)-(13). We find that the ratios of 
two successive binding energies obey the universal scal- 
ing law in Eq. ([!]) expected for the Efimov effect. This 
result establishes that the observed bound states of three 
magnons arc indeed the Efimov states. Note that it is in 
general possible that there are also non-Efimov three- 
magnon bound states in other channels, which may have 
lower binding energies. 

In the case of s = 1/2, the two lowest binding ener- 



gies near the two-magnon resonance (111 are shown in 



FIG. 3. (Color online) Lowest binding energy of three 
magnons for s — 1 with J z = J (upper panel) or J z = — J 
(lower panel) near the two-magnon resonance ( |13[ ). See the 
caption of Fig. [2] for other details. 



of Efimov states is completely characterized by the s- 
wave scattering length a s and the so-called Efimov pa- 
rameter k* [I]: 



EL -> -A 



-2n 



*-F[ A " 



(n 



oo 



(18) 



Here 1/m = 2s J a 2 is an inverse mass of constituent 
particles and F{x) is the universal function defined in 
a range —0.663293 < x < 14.1314 and normalized as 
F(0) = 1 [4'.. By matching Eq. ( 18 ) with the binding en- 
ergy of the second excited state at the resonance a s — >• oo; 
E 2 /J ~ -8.08 x 10" 7 , we obtain K*a ~ 0.463. The re- 
sulting universal curves for n — and n — 1 are also plot- 



Fig. [2] The universal theory predicts that the spectrum 



ted in Fig. [2] by using a s obtained in Eq. (10). We find 
an excellent agreement of the binding energy of the first 
excited state (n = 1) with the universal theory, which 
leaves no doubt that this three-magnon bound state is 
the universal Efimov state. On the other hand, the bind- 
ing energy of the ground state (n = 0) deviates from the 
universal theory away from the two-magnon resonance 
where non-universal corrections become non-negligible. 

Similarly, in the case of s = 1 with J z = J (or — J), we 
obtain ~ 0.173 (or 0.0478) by matching Eq. Q with 
the binding energy of the first excited state at the reso- 
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FIG. 4. (Color online) Critical exchange anisotropy for the 
two-magnon resonance a„ — > oo for s = 1/2 (solid) and s = 1 
(dashed curve) with D = in the presence of the spatial 
anisotropy 7. An s-wave bound state at K = appears in 
the upper side of the critical curve. 



FIG. 5. (Color online) Critical spatial anisotropy for the two- 
magnon resonance a s — > 00 for s = 1/2 in the presence of the 
frustration J z / J 2 * . An s-wave bound state at K — (0, 0, 2ko) 
appears in the lower side of the critical curve. The critical 
spatial anisotropy for s = 1 with D = is at most 7 < 
3.61 x 10~ 5 and thus invisible here. 



nance a s -> 00; E x /J ~ -1.16 x 1CT 4 (or -8.88 x 1CT 6 ). 
The resulting universal curve for n = is plotted in Fig. [3] 
as well as the lowest binding energy near the two-magnon 



resonance (13 1. We find a reasonable agreement between 



them which confirms that this three-magnon bound state 
is consistent with the universal Efimov state. 



REMARKS ON EXPERIMENTAL REALIZATION 

In summary, we showed that magnons in quantum 
Heisenberg magnets exhibit the Efimov effect by tun- 
ing an easy-axis exchange or single-ion anisotropy. The 
single-ion anisotropy D can be changed significantly 
in organic magnets by choosing different ligands [40] . 
Therefore, it is possible to find a compound with spin 
s > 1 whose D/J ratio is already close to the critical 
value. For example, there is an s = 1 ferromagnetic 
compound based on Ni with D/J ~ 3.0(5) [41], which 

Further- 



is not far from the critical value in Eq. (13a 



more, the exchange coupling J can be tuned by applying 
pressure [42] to bring the system near the two-magnon 
resonance to realize the Efimov effect. 

On the other hand, the single-ion anisotropy plays no 
role for s — 1/2. Therefore, the two-magnon resonance 
needs to be induced by the exchange anisotropy. How- 
ever, the exchange anisotropy is in general difficult to 



tune and its critical value in Eq. (Ill is somewhat large 



Here we point out that the critical exchange anisotropy 
J z / J can be reduced significantly if a magnet is spatially 
anisotropic. The spatial anisotropy of the exchange cou- 
plings can be taken into account simply by replacing J 
and J z in the Hamiltonian ^ with J e and JJ, respec- 
tively. We assume a uniaxial anisotropy J = J x = J^ 
and J z = jf = Jf with ratio 7 = J / J z = J z /J*. The 
critical exchange anisotropy obtained by solving Eq. Q 
is shown in Fig. [4] and we find that it reduces signifi- 



cantly toward an isotropic point J z / J — > 1 in a quasi 
one-dimensional (7 — > 0) or two-dimensional (7 — > 00) 
limit. Because magnets with strong spatial anisotropics 
are very common [43H49] , there is hope to find an s = 1/2 
ferromagnetic compound whose J z /J ratio is already 
close to the critical value. Once such a compound is iden- 
tified, the spatial anisotropy 7 can be tuned by applying 
pressure in one or two direction(s) to bring the system 
near the two-magnon resonance to realize the Efimov ef- 
fect. 

An alternative approach to induce the two-magnon 
resonance in s = 1/2 magnets even without the ex- 
change anisotropy J z — J is to introduce a frustra- 
tion. The simplest example is given by quasi onc- 
dimensional spin chains with nearest-neighbor ferromag- 
netic and next-nearest-neighbor antifcrromagnetic cou- 
plings; i?intra = ~ ^2 r (J Z Sr ' S r+ £ ~ J 2z S r ■ S r+ 2z), 

which are realized by Rb2Cu2Mo30i2 [50], LiCuV04 [51"] . 
and Li 2 Cu0 2 [55]. When J z /J 2z < 4, the single-magnon 
dispersion develops a minimum at nonzero momentum 
ko = arccos( J z /4J 2z ), and accordingly, two magnons 
form a bound state with 2fco [S"2Tt58| . Then, by in- 
creasing an interchain coupling; -ffi n t cr = — J2r J (S r ■ 
Sr+x + S r ■ S r +y), the two-magnon bound state disap- 
pears into the threshold, which leads to the divergence of 
an s-wave scattering length. The corresponding critical 
spatial anisotropy 7 = J/J z is shown in Fig. [5] Be- 
cause the frustration ratio J z /J 2z is highly tunable un- 
der pressure due to the strong dependence of ferromag- 
netic couplings on the cation-anion-cation angle of the 
superexchange path jSHj, the Efimov effect can be real- 
ized in s — 1/2 frustrated magnets without the exchange 
anisotropy. The same approach can be used also for quasi 
two-dimensional frustrated magnets [60] , 

Finally, the emergent Efimov states of magnons can 
be observed with previously-used spectroscopic measure- 
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merits [5H30). such as far- infrared absorption, inelastic 
neutron scattering, and electron spin resonance, provided 
that the U(l) symmetry of the Hamiltonian to conserve 
J2 r Sr ^ weakly violated H3 EH]. Our study thus 
opened up new avenues for universal few-body physics in 
solid state systems. 

This work was supported by a LANL Oppenheimer 
Fellowship and the U. S. DOE contract No. DE-AC52- 
06NA25396 through the LDRD program. 
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